Effect of assisted hopping on the formation of local moments in magnetic impurities 

and quantum dots. 
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Assisted liopping effects in magnetic impurities and quantum dots are analyzed. The magnitude 
of the assisted hopping term in a quantum dot in the limit of large level spacing is comparable to 
other corrections induced by the electron-electron interactions. Assisted hopping leads to differences 
between conductance peaks associated to the same level, and, when the effect is sufficiently strong, 
to local pairing correlations. 
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I. INTRODUCTION 

Electron-electron interactions lead to the formation of 
local moments andr-to the Kondo effect in magnetic im- 
purities in meljalsaa, and in quantum dots attached to 
metallic leadsBofl. This effect is the direct consequence 
of the Coulomb repulsion between electrons which occupy 
the same quantum level of the impurity or quantum dot. 

If we assume that the main physical features of the 
impurity (or the dot) are governed by a single quantum 
level, the leading correction to the intralevel Coulomb 
interaction is an assisted hopping term. This interaction 
induces a dependence of the coupling of the impurity (or 
the dot) and its environment on the occupancy of the 
leveE. 

In the following, we estimate the magnitude of such a 
term for models of large atoms and quantum dots, and 
analyze its effects on the formation of a local moment 
on the impurity, or quantum dot. For bulk systems, the 
inclusion of an assisted hopping term in the electronic 
hamiltonian favors the existence of pairing correlationsQ. 
In the case of an impurity, this tendency towards local 
pairing quenches the local moment, and, for quantum 



dots, it can lead to asymmetries in the conductance of 
peaks associated to the same level, and even to an en- 
hancement of the conductances at low temperatures. 

The model is described in the following section IIA. 
Simple estimates of the magnitude of the assisted hop- 
ping term for quantum dots are discussed in section IIB. 
A mean field analysis is given in sections III A and IIIB. 
Extensions of the mean field approach are presented in 
section IIIC. Section IV contains the main conclusions. 



II. THE MODEL 



The hamiltonian. 



We analyze a single quantum level, associated to the 
creation operator c?|, where s =t, i is the spin. This level 
is coupled to a continuum of non interacting electrons, 
described by operators c\^. For simplicity, we assume 
that a single-channel in the environment interacts with 
the impurityEI. The hamiltonian is: 



J 



o 
o 



(1) 



where V is the volume of the system, and: 



ndi 
nd 



nd-\ + ndi 



(2) 



The density of states of the conduction band at the Fermi 
level is N{ep). We assume that VN{eF), /WN{ep) < 1. 
Without loss of generality, we will set = 0, and define 
No = N{eF). 



When AV = we recover the standard model pro- 
posed by Anderson forJhe study of the formation of lo- 
cal moments in metalso. This model is determined by 
three parameters, the Coulomb repulsion U , the position 
of the level e°, and its width, F = V'^N^- The model has 
electron-hole symmetry around = e[] — C/ /2. 

The assisted hopping term in eq.(|l]) is the next leading 
interaction in terms of d operators which can be added 
to the Anderson hamiltonian, assuming the same trun- 
cated Hilbert space. As discussed later, it can be derived 
from intradot interactions only. The hamiltonian (|^) in- 
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torduces an additional dimensionless parameter, IS.V/'SL. 
For small atoms, ISV/V is a number of order unityO. 
The inclusion of this term in the hamiltonian breaks the 
electron-hole symmetry of the initial Anderson hamilto- 
nian. 

III. CALCULATION OF AV. 

A. Definition of AV . 

We assume that the assisted hopping term in eq.(||) 
arises from interactions which involve electronic levels 
within the dot only. When the hamiltonian is trun- 
cated to the Hilbert space defined by a single state within 
the dot, and the leads, the electron-electron interactions 
within the dot can induce, among others, a term of this 
type. 

Our initial hamiltonian is: 

'H — Hdot + "Hlead + "Htunn 

T-idot = ^ Cicjci + X! ^ijkiclcjclci 

i ijkl 

OG 

Ti-lead = ^ X! ^n^"+l + 
n=0 

'Htunri = ^ Vic\c^) + h.C. (3) 

i 

We have assumed that the lead contains a single chan- 
nel. The index i labels all quantum numbers of the states 
within the dot, including spin. We are interested in tran- 
sitions where the charge state of the dot changes from 
iV to iV + 1, and from iV + 1 to iV + 2. In the absence 
of interactions, they correspond to the filling of a par- 
ticular level, which we denote i. When the charge state 
of the dot is N, all levels j such that j < i are occup- 
pied. We use as basis for the electronic states in the dot 
the Hartree-Fock wavefunctions defined when for charge 
state N. Rearrangements of the electronic levels imply 
that the corresponding waveiiinctions for charge states 
+ 1 and A^ + 2 are differentEI. Typically, the tunneling 
between the dot and the lead takes place in a region of 
atomic size at some point of the surface of the dot, Fq. 
Then, Vi cx ^^(fo), where 5'i(f) is the wavefunction of 
state i. 

We now assume that the broadening of the levels due 
to the coupling to the leads, T = (V^Nq) is smaller than 
the mean level spacing, A. Then, the effective tunneling 
can be estimated from the ground state wavefunctions of 
the dot for different charge states: 

V:'^J^''+' = {N + l\Y,V,c]\N) (4) 
j 

and the difference between adding two electrons to state 
i, and adding one is: 

AV = - yjr^+i (5) 




|N+l)^ |N+l)^ |N+l>^ 

FIG. 1: Sketch of the zeroth order wavefunctions, and correc- 
tions with one and two electron-hole pairs for charge states 
and A'^ -f 1. See text for details. 



B. Perturbative analysis. 

For large dots, the interactionJ-emas decrease with the 
size and conductance of the dotEHltirllj, so that a pertur- 
bative calculation of AV is possible. 

As mentioned above, we use the electronic basis which 
diagonalize the Hartree-Fock approximation to Tidot in 
charge state A^, which we denote Ti^p- We can define 
"Hint = 'Hdot — 'H^p, and the wavefunctions: 

l^>o - \N)hf 
\N + l)o = 4^\N)hf 
\N + 2)o - 4^cI\N)hf (6) 

where |A^)ifF is the ground state of H^p (note that |A-t- 
l)o and I A^-|-2)o are not the Hartree-Fock approximations 
to the wavefunction of the dot with charge A^ + 1 and 
N + 2). 

The interaction hamiltonian, 7ii„t, when acting on 
|A^)o induces two electron-hole pairs. The occupation 
of level i in wavefunctions |A^-|- l)o and \N + 2)o, and the 
induced interaction between this electron and the rest im- 
plies that the basis which diagonalizes H^p is not longer 
optimal. Then, Timt acting on these wavefunctions leads 
to excited states with one and two electron-hole pairs. A 
sketch of the corrections to the wavefunctions is shown 
in Fig.fll. 

Matrix elements which involve one and two electron- 
hole pairs have a different dependence on the size of the 
dot (see below). We will calculate Ay to second order in 
the matrix elements associated to one electron-hole pair, 
'Hint, I, and to first order in the matrix elements involving 
two pairs, 'Hint,2- We use eqs. (^ and (||), and we need 
to calculate the wavefunctions |A), |A -f 1) and |A^ + 2). 
Generically, we can write: 
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^ f 1 - i y I('^iI^'"m|0)IM ( ^ y |ni)(ni|H.„M|0) ^ y \n2){n2\nint,2\0) ^ 

V- |n2)(n2|H,„t,i|^i)("i|^^nt,i|0) V- |ni)(ni|H,„t,i|0)(0|Hznt,i|0) \ , . 



where |0), |ni) and \n2) describe states with zero, one and 
two electron-hole pairs. 

The operator J^i ^i^l eq.(^) can, at most, change by 
one the number of electron-hole pairs in the wavefunc- 
tion. Hence, there are no contributions to eq.(^) which 
are of first order in Ti.int.2- In addition, the terms in the 
second line in eq.(0) are of second order in Hint,i and 
lead to wavefunctions with, at least, one electron-hole 
pair. Their contributions are of third order in Timt,!, 



and need not be considered. As the basis set used gives 
the Hartree-Fock solution for charge state N, we have 
that: 

n,ntA\N)HF^O (8) 

which implies that the corrections to V^J*^^^ are of 
second order in TimtA- Finally, adding all contributions, 
we find: 



AV: 



E 

k>i 



k>i I 2 

oy. V 



k>i 



Vjhjkhik 



k>i 

E 



Vkhjkhik 



- efc)(ei - efe) ^ (e^ - - ek) 



(9) 



where: 



hjk = (j|Hi„t,i|^) 



(10) 



and we assume that all levels are doubly degenerate. 

Eq.(^ implies that electron-electron interactions 
within the dot give rise to a first order correction to the 
hopping, whose sign depends on the nature of the in- 
teraction, and a second order correction which tends to 
suppress the hopping. In the limit where the level spacing 
within the dot is negligible with respect to the tempera- 
ture and charging energy, these terms change into the non 
equilibrium correction.s to the effective tunneling density 
of states analyzed int^EJ. Note that in in addition to the 
correction to the hopping given in eq.(^), the polariza- 
tion of the occuppied levels, j < i, gives a correction to 
V which does not depend on the number of electrons in 
state i. 



C. Spherical dot. 

We now apply the previous analysis to the case of a 
spherical dot of radius R and iV 1 electrons. For 
simplicity, we assume that the positive charge needed to 
stabilize the system is uniformly distributed within the 
radius R. The level spacing is: 



mR 



(11) 



where kp — (97r/4)^/^i?~^ is the Fermi wavevector, and 
m is the mass of the electrons. For sufficently large elec- 
tronic densities, kp 3> (rne'^)/h'^, the Coulomb repulsion 
between electrons can be treated using the Fermi- Thomas 
approximation. Within this approximation, Timt.i is 
given by the Fermi-Thomas potential induced by the ad- 
dition of a unit of of charge to the dot. |— .For a three 
dimensional spherical dot, this potential ialSl: 



VpTir) 



^2^-kFT{R-r) 



kpriR-ry 



(12) 



where 



is the radial coordinate, and kpT 



(Ae^ mk f) / {t^Tt- ) is the Fermi-Thomas wavelength. In 
addition to this term, we have to include a short range 
interaction which leads to the excitation of two electron- 
hole pairs when expanding around the Hertree-Fock so- 
lution, 7iint,2 in eq.(^. The matrix elements of TLmt^i 
decay as R~^ for kpxR S> 1, while those of 71^,2 de- 
cay as R~^, justifying the separation of scales used in 
obtaining eq.(||). 

The potential VpT is localized within a shell of size kp}p 
around the edges of the dot. We assume that the matrix 
elements {k\VpT\j) are roughly constant if Ak <^ kpT 



where e,- 



hvpAk, and zero otherwise. As the 



level spacing is given in eq. ( |11[ ) , the number of levels for 
which the matrix elements of VpT are not negligible is 
kpxR- Within this energy window, we can also assume 
that |^'fe(fo)| ~ l^j(?o)|j so that the the correction to 
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AF, using eq.(^, is the sum of kpxR terms of the same 
sign. Finally, Vpx conserves angular momentum. The 
maximum angular momentum in a sphere of radius R is 
approximately kpR. 
Defining: 



Vft = (»|VftK) 



kprpR'^ mkpR^ 



(13) 



we find that the first (linear) and second (quadratic) con- 
tributions to ISV are: 



A^i cx y-^log(fcFTi?) 



A 

IW2 oc VikpK) 



2 ^ FT 
A2 



logikprR) 



(14) 



where the {kpR)^ factor to AV2 arises from the number 
of angular momentum channels. We have obtained that 
the two corrections are of similar magnitude, and decay 
as (kpR)'^- 



As in the case considered in the previous subsection, 
the leading term in an expansion in arise from the 
Fermi- Thomas potential, eq.(12), upon the addition of 
electrons. This potential is determined solely by the ge- 
ometry of the dot and electrostatic constraints, and is 
independent of the details of the dynamics of the elec- 
trons. 

We use the_matrix elements of Vft calculated for a dif- 
fusive dot iiiE£l. As in eq.(p^), we find two contributions, 
AVi and AV2. The average over disorder of AVi is zero, 
with mean fluctuations: 
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(15) 



where ci is a dimensionless constant of order unity. The 
log{g) correction is due to the summation over states such 
that A < e,- - efe < Et- 



To next order in g ^, we find: 

AV2 cc V- log(g) 
9 



(16) 



D. Diffusive dot. 

The low energy levels and wavefunctions of quantum 
dots in the diffusive regime are well described using ran- 
dom matrix theory. We consider a dot of average ra- 
dius R, mean free path I, and mean level separation 
A. In addition, we assume that the electron density 
is large, < 1, where rj^ = {4/9TTy/^{h'^kF)/{me'^). 
The electron-electron interaction can be expanded in 
powers of the inverse-jpriductance, g = Et/A, where 
Et = {n^kFl)/{mR^ 



For dots of similar radii, AT^ is larger in the diffusive 
regime than in the regular case. For chaotic, ballis- 
tic dots, it seems reasonable to replace g by kpR in 
eqs.(|l|Jll). 

IV. MEAN FIELD ANALYSIS. 

A. Magnetic solutions. 

We now consider the hamiltonian in eq. (0). We first 
analyze possible magnetic solutions. The mean field 
hamiltonian is: 



Hmf = y^gfc4^^Cfc,^ + egnrf + U {{nd\) Udi + {ndj) fid]) + Vy^^cl^ds + h.. 



AV 



(17) 



r 



where cq is defined in eq. ^ . The sites defined by Ci , i 7^ 
can be integrated out. leading to a 2 x 2 matrix equation 
for the components of the Green's function for each spin 
index projected on the cq and d sites: 



-V + AVns 



-V + AVus ~ ej - Ung + AVg^ 



where s is the spin index, and we have defined: 



[nds 



(18) 



(19) 



For the nonmagnetic solution, the different components 
of the density of states are (omitting spin indices): 



ImGdd(w) 
ImGoo(w) 
ImGodiuj) 



{uj ~ edf + r2 
N^{LO-edf 



where: 



F 

Erf 



{lu - edf + r2 
jVo(^-£rf)(F- AT/np) 
(c. - erf)2 + r2 



Na[V~AVn^f 
eO + C/no - AVg^ 



(20) 



(21) 
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The self consistent relations required in eq.(|T^, and the 
Green's functions in eq.(po|) imply that: 



no 



1 1 /Ed 

arctan — 

2 TT vr 



27r 



log[NSiel + r')] (22) 



For the symmetric Anderson model, we have Cd = 0, 
and hq — 1/2. The equations for the magnetic solution, 
in the limit of small magnetization, are analyzed in the 
Appendix. 

Numerical results of the occupancies of the different 
spin states as function of U, V, Ay and eo are shown in 
Fig.g. 



FIG. 2: Number of electrons per spin, as function of U. Dot- 
ted curves: AV = 0. Broken Curves: AV = 0.15. In all 
cases, V = 0.2, in units of N^J- ■ Top: + U/2 = 0. Center: 
+ U/2 = -0.05. Bottom: + U/2 = +0.05. 



B. Superconducting solutions. 



Alternatively, we can use the BCS decoupling and 
write a mean field hamiltonian: 



J 



Ti-Bcs = y^te^ .jCn+i.s + edUd + V^y^ + c^Icq^.s'I + U (did] ) d^d 

n,s s 



I 



h.c. 



(23) 



This hamiltonian couples up spin electrons with down on the electron and hole states at sites cq and d can be 
spin holes, and viceversa. The Green's function projected written in terms of two 4x4 matrices: 

I 



-V 


\ AVn 



-f-d 
AVn 
^Un + AVg 

r 





AVn 
V 



AVn \ 
-Un + AVg 
V 

+ ed 



(24) 



where the two upper (lower) rows correspond to electrons 
(holes), and we have defined: 



Corrections beyond the mean field 
approximation. 



n = {d-\di) 

9 = {d-\ci+dic^) 



(25) 



The selfconsistency condition implicit in eq.(|24|) can be 
linearized, as outlined in the Appendix. 



The impurity model studied here can be considered as 
a (0+1) dimensional model, where fluctuations in time 
need to be considered. In principle, we can integrate 
out the fermion in the leads, and obtain an effective four 
state model. The four states correspond to the four occu- 
pancies of the single state at the impurity, as defined in 
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u 



FIG. 3: BCS order parameters, eq.(||), for V = 0.2, = 
0.15, and ed ~ U /2 — Q va units of TV^ef)^^. The transition 
is weakly first order. Broken line: occupancy of the localized 
level (right scale). 



that there are two physical regimes: i) Near rid = 1 and 
for large values of U /T , a local magnetic moment is in- 
duced, and, ii) away from Ud — \, and for moderate 
values of CZ/F, pairing correlations develop near the dot 
or impurity. 

We identify the local moment regime studied earlier 
with the Kondo effect. Note, however, that the relation 
between the Kondo temperature and the parameters in 
the hamiltonian will be changed by the assisted hopping 
term. 

Pairing fluctuations imply that charge states which dif- 
fer by two are strongly mixed in the ground state. This 
can be checked by the exact diagonalization of a block 
containing the impurity and a few sites in the leads. The 
assisted hopping term leads to an effective negative Uef / 
model, where C4// = Ej^ -\- — 2Enj^i, for certain 

values of N . These relatively small blocks, when attached 
to the rest of the leads will act as local defects with a ten- 
dency towards pairing. Hence, we conclude that the BCS 
solution identified in the mean field approximation corre- 
sponds to a regime described by the negative U Anderson 
model. 
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FIG. 4: Mean field phase diagram of the hamiltonian in eq.(|l|). 
SC: superconducting solution. The thick line is a first order 
transition (see Fig.[H]). The jump in is comparable to the 
thickness of the line. LM: local moment. For comparison, the 
results for = are also plotted. 

eq.(0). These states have long range interactiopjiij time, 
as in related fermionic and bosonic problemsll3ll3. The 
fluctuations in time will tend to restore the symmetries 
broken at the mean field level. We will use the mean field 
results as a guide to infer which type of universality class 
will determine the low temperature properties. The fact 
that the reduced Hilbert space associated to the impurity 
includes only four states greatly restricts the possible low 
energy fixed points. 

The analysis in the two previous subsections suggests 



V. CONCLUSIONS. 

We have analyzed the influence of assisted hopping in 
correlated impurities embedded in a metal, or quantum 
dots coupled to leads. This type of interaction is the 
simplest term which can be defined in the restricted basis 
given by a single quantum level within the dot. Related 
effects, in the opposite limit where the states withiaJie 
dot can be treated as a continuum were studied inOEj. 

We have estimated, in section II, the magnitude of this 
term for simple models of regular and diffusive dots. In 
the latter case, we find that the assisted hopping term has 
a similar dependence on the conductance of t]a«.|dO|t-as-jtiie 
interaction corrections to the peak spacingl30't3't3c3. 
The assisted hopping term is directly related to the 
changes of the wavefunctions within the dot upon the 
addition of a single electron. This interaction will be 
larger than estimated in section II in one and two di- 
mensional geometries as the Fermi- Thomas potential is 
more extended, and for low electronic densityO. Note 
also that the orthogonality catastrophe is expected,to be 
enhanced in disordered two dimensional systcmsEll. We 
have only taken into account the interactions within the 
dot. In devices with more than one dot, the interactions 
between electrons in different dots will also enhance the 
effects reproted here. Assisted hopping can contribute to 
change the peak height distribution with.respect to that 
predicted by Random Matrix Theoryn3c3. 

Our results suggest that assisted hopping can lead 
to local pairing correlations for reasonable values of 
the parameters. This implies enhanced conductivity 
through the dot or impurity, significant deviations from 
the electron- hole symmetry implicit in the Kondo effect, 
and measurable differences in the conductance through 
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the two peaks associated to the same quantum level. The gn « go ± 5g (26) 

existence of these effects does not seem incompatible with 
present experimental evidencec^. A detailed analysis of 
assisted hopping effects quantum dots, in the limit of 
large level spacing, deserves further investigation. 
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The Hartree Fock approximation introduces four varia- 
tional parameters, e^i, Cdj., Tj and F^. We also make the 
expansion: 



VII. APPENDIX. €dii « eQ±Se 

Tu « To±5r (27) 

A. Linearization of the magnetic Hartree-Fock 
equations. 

We first analyze the magnetic solutions when the mag- 
netization is small. Then, we can write: 



Hdii ~ no±Sn The consistency between eqs.(E6|) and eqs.(27) leads to 



6n 

Sg 

Se 
5T 



ToSe 



eoST 



NqAVSu 



log 



27r 

-USn + AVSg 
2NoVAV6n 



Eqs. (psl ) have a non trivial solution if: 



Det 



To 





N„{V-n„AV)eo 



NojV- AVn„){eoSe + ro6r) 



u 
1 

-2AroVAV 
NqAV , / w\ 







1 

No{V-noAV)r„ 



~AV 




1 



(28) 



(29) 



For AV = 0, eq. 



|) reduces to: 



For eo = (symmetric case), and AV ^ 0, we find: 



C/ = ^Fo 



21/ 



V - n^AV 



log 



(30) 



(31) 



This equation defines the critical value of U at which a 
solution with a non zero magnetization appears. 



B. Linearization of the BCS equations. 



We define Gdd{^^) and Gdcgi^) as the anomalous 
Green's functions involving electrons and holes induced 
by the BCS coupling in eq.( P^). They can be obtained 
from the implicit equation (p4|) . The selfconsistency 
equations are: 



n = — 

TT 



lmGdd{uj)duj 
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5 = - 



(32) 



We can invert eq. (g4|) and linearize with respect to n and 
g, to obtain: 



J 



ImGrfd(w) 
ImGrfco (^) 



2yAyAro(w2 -el- r^)n + 2rw(C/n - AVg) 
2r2AT/iVown + _ g2 _ p2) [7VgAV^(w - ed)n - NoV{Un - AVg)] 



[(c. - e,)2 + r2] [(a; + 6^)2 + r2] 



(33) 



We can integrate these expressions, so that: 

Un - AVg 



arctan 



m 
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^^^{log[A^o^(.^ + r^)] + arctan f!^n 



T J ) 



(34) 



Finally, we obtain: 



arctan f — ^ 

ntd V r / 



. arctan (^) log [iV„2(e^ + n] 



< 



(35) 



where we have dropped the arctan term in the second 
equation in (pj 
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